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Abstract 

(N 

^ • In May 2002, the BL Lacertae object (BL Lac), lES 1959+650, was detected in 

2 ■ the Very High Energy (VHE) regime by the VERITAS collaboration using the Whipple 

7-ray observatory. During the following two months, the VERITAS collaboration ob- 
^ ■ served episodes of flaring; at times, the nightly average VHE signal was three times as 

O ■ bright as the Crab Nebula, the standard candle in VHE astronomy. lES 1959+650 was 

also monitored in the optical and X-ray bands. The X-ray light curve showed periodic 
behavior and correlation with the VHE and the optical bands. 



O 
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^1 1. Introduction 

a 

The BL Lac object, lES 1959+650, redshift 0.047, is an active galaxy pow- 
^ ■ ered by a super massive black hole. It was discovered in 1993 by comparing the 

■ X-ray/radio/optical fluxes of objects in the Einstein IPC Slew Survey (Schachter et. 

al. 1993). 1ES1959+650 was first detected at VHE energies by the Utah Seven Tele- 
scope Array collaboration (Nishiyama et. al. 1999) in mid- 1998. It was seen in a 
flaring state in the VHE energy regime, this Spring for the first time by the Whipple 
Observatory (Dowdall et. al. 2002). Following the initial detection, the object was in- 
tensely monitored by the VERITAS collaboration during the next two months (Holder 
et. al. 2002a and Holder et. al. 2002b). Overlapping with the VHE observations were 
R-band photometry, optical spectroscopy, and X-ray observations. The study of the 
cross-correlation of these energy bands is an important tool in restricting the possible 
emission and absorption processes of VHE photons. 
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Fig. 1. Light curve of lES 1959+650 in the VHE, X-ray, and optical energy bands averaged 
over 24 hours. 



2. Observations and Data Reduction 



2.1. VHE observations 

The Whipple Observatory 10 m reflector images Cerenkov light flashes from 
air showers initiated by VHE photons in the atmosphere. The telescope is located 
on Mt. Hopkins near Tucson, Arizona. It consists of 238 mirror facets arranged on 
a spherical support structure of Davis-Cotton design, with focal length 7.3 m. The 
camera consists of 490 ultraviolet- sensitive photomultiplier tubes, giving a 3.8 ° field 
of view. Only the inner PMTs covering 2.6 ° are used in this analysis. Observations on 
lES 1959-1-650 were taken from 16 May 2002 through 8 July 2002 (UT) for a total of 36 
hrs. ON source and 6.4 hrs. OFF source. The detection significance for ON/OFF paired 
observations is 8 a with an average photon rate of 1.5 ± 0.2 7/min. Standard analysis 
of the data was done by moment fitting (Hillas 1985) of the events followed by a set 
of cuts optimized on Crab Nebula observations to optimize the signal-to-noise ratio. 
The light curve Fig[l] was calculated from ON source data, i.e. tracking analysis, with 
background estimated from OFF source data. The rate was then corrected empirically 
for elevation and throughput (LeBohec and Holder 2002) based on observations of the 
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Fig. 2. Top: ASM X-ray light curve, 1 a errors indicated by dots. Bottom: Fourier transform 
of light curve. 

Crab Nebula, the standard candle at VHE energies, taken from February to April 2002. 
Strictly, this is only valid if the sources have the same spectral index 7; 7crafe = 2.5. 

2.2. X-ray observations 

X-ray measurements were taken by the ASM detector on board the RXTE satel- 
lite. The ASM quick-look results were provided by the ASM/RXTE team. They cal- 
culated the 24-hour average flux from five to ten 90-second dwells. The ASM detects 
photons in the energy range from 1.5 keV to 12 keV, with peak sensitivity at 5 keV. 

2.3. Optical observations 

The R-band (650 nm) photometry was taken with the 70 cm telescope at the 
Abastumani Observatory in Georgia from 19 May to 11 June 2002. The frames were 
reduced using DAOPHOT II. The apparent magnitude of IBS 1959-1-650 was deter- 
mined by comparison with 3 standard stars in the field of view (Villata et. al.l998). 
Optical spectra were taken using the newly refurbished 6.5 m telescope at the MMT. 
The spectra were taken while 1 ESI 959-1-650 was active at VHE energies and one week 
later when the flux was not detectable. The first of the two spectra was measured on 
5 June 2002 with a 300 lines/mm grating and 45 min exposure, the second spectrum 
was taken on 13 June 2002 with a 500 lines/mm grating and 10 min exposure. Not all 
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cosmic ray contaminants were removed from the spectra. 
3. Light curve correlations 

3. 1 . Quasi-periodic X-ray flaring 

The X-ray light curve, Fig [2], shows a series of flares spaced roughly 40 days 
apart. Also shown is the associated Fourier transform, Fig[2]. It shows that there 
is weak evidence for a 40 day period of X-ray flaring. A similar pattern of flaring 
activity was possibly found for Mkn 501 (Kranich et. al. 1999). 

3.2. VHE and X-ray correlation 

The X-ray - VHE linear correlation coefficient (LCC) as a function of the time 
delay between the bands is shown in Fig[3]. Simultaneous keV-VHE measurements are 
correlated with a LCC = 0.45. Here, simultaneous means ±12 hours. To determine the 
significance of this correlation, the keV and VHE flux measurements were randomized 
100 times with respect to one another and the LCC was calculated for each random set. 
This resulted in a mean LCC = 0. 1 ± 0. 1 , so we deduce that the correlation between X- 
ray and VHE flux is significant at the level of 3.5 a. Assuming the LCC distribution to 
be Gaussian, the probability of at least one 3.5 a deviation in 61 independent samples 
is 0.03. This is an approximation as neither the LCC distribution is Gaussian nor the 61 
time intervals are completely independent samples. On the other hand, when the X-ray 
flux is time delayed by 25 days following the VHE flux, we find an anti-correlation of 
the X-ray and VHE flux. The LCC factor is 0.42, the probability for this to happen 
by chance is 0.15. So, there is weak evidence for both a 40 day periodicity of the X- 
ray flux and simultaneous correlation with VHE photon emission around the time of 
observation. A simultaneous correlation of X-ray and VHE flux has also been found 
during flares of Mrk421 in 1995 and Mrk 501 in 1997 (Buckley 1999). 

3.3. Optical and X-ray correlation 

In the optical, lES 1959-1-650 was very dim compared to the light curves from 
1996-97 (ViUata et. al. 2000) and 1997-99 (Kurtanidze et. al. 1999). During those 
times the magnitude varied from 15.05 to 14.55. Thus, during our VHE detection the 
source was optically fainter and not as variable as during previous observations. The 
X-ray and optical (R-band) emission appear to be correlated, although with the latter 
delayed or lagging by 5 days. The LCC is 0.78, 5 a away from a Gaussian random 
distribution. The probability for this to occur randomly in 16 independent samples is 
9 X 10"^. Delayed correlation between R-band and X-ray emission has also been seen 
in a flare of Mrk421 in 1996 (Buckley 1999). 
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Fig. 3. Correlation of X-ray and VHE rate. 



4. Optical spectra 

Typical of a BL Lac object, the optical spectrum of lES 1959+650 was nearly 
featureless. The spectra were neither flux calibrated nor corrected for atmospheric 
extinction. The relative change in spectral features was calculated by first normalizing 
to the smoothed continuum and then subtracting the two normalized spectra from one 
another. The smoothed continuum was calculated by first applying a 3 pixel wide 
boxcar average to the 300 lines/mm spectrum and a 5 pixel average to the 500 lines/mm 
spectrum to minimize any effects the different gratings had. Then, the continuum was 
calculated as a 100 pixel (200A) boxcar average. The relative change in the optical 
spectrum between states of high VHE emission and low VHE emission was measured 
after applying a 3-pixel box smoothing to reduce the sensitivity to noise and possible 
mismatched wavelength calibration. The intensity of the following lines showed a 
decrease from 5 June to 13 June: if„ (gfwhm = 9 A) 3%, Mg (gfwhm = 9 A) 2.5% 
, and CaFe (gfwhm = 11 A) 2 %. Unfortunately, the precision of the measurement 
for the Ha is uncertain because it is located in a region of strong night-sky absorption 
lines. 
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Fig. 4. Left: The dependence of the R-band - VHE linear correlation coefficient on the time 
delay between the bands. Right: Correlation with a 5-day delayed optical magnitude. 

5. Summary 

The BL Lac object, lES 1959+650, was detected in Spring 2002 with the Whip- 
ple 10m Cherenkov Imaging Telescope. The 24-hour average VHE photon flux was 
measured to vary by as much as 3 Crab units. The VHE and X-ray flux showed a slight 
simultaneous correlation. There is evidence for a roughly 40 day period in the X-ray 
light curve. A strong correlation between the 5-day delayed optical R-band brightness 
and the X-ray flux was also observed. The R-band brightness was fainter than dur- 
ing previous observations and optical spectra were without any prominent features or 
variability. 
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